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ABSTRACT 

A possible accretion model associated with the ionization instability of quasar 
disks is proposed to address the growth of the central black hole harbored in the 
host galaxy. The evolution of quasars in cosmic time is assumed to change from a 
highly active state to a quiescent state triggered by the S-shaped ionization insta¬ 
bility of the quasar accretion disk. For a given external mass transfer rate 
supplied by the quasar host galaxy, ionization instability can modify accretion 
rate in the disk and separates the accretion flows of the disk into three different 
phases, like a S-shape. We suggest that the bright quasars observed today are 
those quasars with disks in the upper branch of S-shaped instability, and the 
faint or ’dormant’ quasars are simply the system in the lower branch. The mid¬ 
dle branch is the transition state which is unstable. We assume the quasar disk 
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evolves according to the advection-dominated inflow-outflow solutions (ADIOS) 
conflguration in the stable lower branch of S-shaped instability, and Eddington 
accretion rate is used to constrain the accretion rate in each phase. The mass 
ratio between black hole and its host galactic bulge is a nature consequence of 
ADIOS. Our model also demonstrates that a seed BH (~ 2 X IO^Mq) similar to 
those found in spiral galaxies today is needed to produce a BH with a final mass 
(~ 2 X 10^Mq). 

Subject headings: accretion: accretion disk - galaxies: active - galaxies: bulge - 
galaxies: nuclei - quasars: general 


1. Introduction 

Nuclear feeding is most important in the evolution of quasars and quasar host galaxies. 
Since quasars were first discovered, it has been suggested to be powered by the accretion of 
gas on to super-massive black holes (SMBHs) at the center of galaxies (Lynden-Bell 1969). 
In the framework of the hierarchical dark-matter (DM) cosmology (Haehnelt et al. 1998), the 
formation and evolution of galaxies and their active nuclei are intimately related (Sanders 
et al. 1996; Boyle et al. 1998; Fabian 1999; Granato et al. 2001; Haehnelt et al. 2000). 
Consequently, the growth and the fuelling of SMBH, and the construction of the galactic 
bulge, are now seen as being potentially related processes (Silk et al. 1998; Fabian 1999; 
Granato et al. 2001). Sub-millimeter photometry of eight X-ray-absorbed active galactic 
nuclei (AGNs)(Page et al. 2001) reveals that active nuclei and the galaxy evolve together, 
with the black hole (BH) accreting matter and the galaxy creating stars. On the other hand, 
recent high-resolution observations of galactic centers have revealed that the estimated mass 
of a central ’’mass dark object” (MDO), which is the nomenclature for SMBH candidate, 
correlates with the mass of a galactic bulge. The BH - bulge mass relation with M^h ~ 
0.001 — 0.006Mbuige (Kormendy et al. 1995; Magorrian et al. 1998; van der Marel 1999; 
Gebhardt et al. 2000a; Ferrarese et al. 2000a) for normal galaxies was shown to extend to 
nearby quasar host galaxies as well (Laor 1998; Nelson et al. 2000; Ferrarese et al. 2001). 
Based on the remarkable dependence of black hole mass on spheroid velocity dispersion 
(Ferrarese et al. 2000a; Merritt et al. 2001; Gebhardt et al. 2000a,b), several theoretical 
models for BH growth have been considered to explain the BH-to-bulge correlations, e.g., a 
hydrodynamic model, including wind regulation (Silk et al. 1998), an inside-out accretion 
model (Adams et al. 2001) and a self-interacting dark matter model (Ostriker 2000) or star 
formation regulated model (Burkert & Silk 2001). As an alternative model for BH growth, 
Umemura (2001) suggested a radiation-hydrodynamical model incorporating the physics of 
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angular momentum transfer. The growth history of the super-massive black holes (BHs) is, 
however, rather uncertain. Two scenarios are involved in modeling these BHs: One is that 
SMBHs can grow out of low-mass ’’seed” BHs through accretion (Rees 1984), and another 
is that SMBHs grow by merging (Barkana et al. 2001). 

It is widely believed that SMBHs built up during the quasar epoch, around 2 : ~ 2 — 3, 
where is the red-shift of quasars. However the observation of optically bright quasars 
shows that if the quasars were continuously radiating, they would grow too large and run 
out of fuel (Haehnelt et al. 1998; Richstone et al. 1998). On the other hand, if the quasar 
disk were unstable, then quasars would be intermittent, and mostly non-radiating, with 
the corresponding long periods of quiescence resolving the old fuelling problem in quasars 
(Shields et al. 1978). This implies that the growth of BHs through accreting gas only happens 
intermittently (Small & Blandford 1992). 

The S-shaped ionization instability hrst found to be responsible for the outbursts of 
cataclysmic variables (CVs) (Meyer et al. 1984; Smak 1984), also operates in binary systems 
containing accreting neutron stars or BHs. The time evolution of the limited cyclic behavior 
is caused by the thermal and viscous instability of the disk in the region of partial ionization 
of the hydrogen and leads to transient behavior of many sources on a time-scale of years 
(King 1995; Cannizzo 1993; Hameury et al. 1998). It is suggested that the similar instability 
may occur in the quasar disk (Lin et al. 1986; Cannizzo 1992; Mineshige 1990; Clarke et al. 
1989; Siemiginowska et al. 1996; Hatziminaoglou et al. 2001), the expected time-scale being 
of thousands to millions of years. Such instability may modify the accretion rate in the 
central part of the flow and thus be responsible for a temporary change of the object status 
from highly active to quiescent on similar time-scales. Burderi et al. (1998) have summarized 
a decade of work on quasar disk instabilities and conclude that all observed active galactic 
nuclei (AGN) are in the outburst phase. For a central black hole with mass of IO^Mq, the 
recurrence time scale of the S-shaped instability is about lO^yr, which the quasar disk spends 
about 10% of the time in outburst phase and 90% in the quiescent phase (Siemiginowska 
et al. 1999). It is thus suggested that the bright quasars we observe today are in a highly 
active state and the faint or ’’dormant” systems are simply those in quiescent states. 

Motivated by the above investigations, we here develop an alternative accretion model 
associated with the S-shaped instability of the quasar disks for the growth of BHs harbored 
in quasar host galaxies. The mass relation between BH and the galactic bulge in which the 
quasar lies is also discussed. Contrary to the former BH growth model through accretion 
(Rees 1984), we address the growth of the BHs using the ionization instability of the quasar 
disk. Eddington accretion rate is used to limit our calculations. The behavior of the disk 
is successively ’’drought”, ’’inundation” and ’’deluge” (Blandford 2002) in the lower, middle 
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and upper branches of the quasar disk, respectively. This paper is organized as follows: 
in section 2, we discuss the accretion features of the quasar disk around a BH in different 
phases, and study the growth of the BH with the evolution of the quasar disk; the BH - 
bulge mass relation is discussed in section 3; hnally, we give our discussion and conclusion 
in section 4. 


2. The growth of BHs associated with the S-shaped ionization instability of 

quasars 

We hrst assume that there is a seed BH in the center of a quasar host galaxy with a 
mass much less than its present day mass, for example, a seed HB could be formed from an 
immediate coalescence of the small-mass BH contained in the merging galaxies (Begelman et 
ah 1980; Cattaneo et ah 1999) or an early formed massive BH with IO^Mq (Umemura et al. 
1993). A massive BH with ~ IO^Mq can also be formed by the collapse of a rotating supper- 
massive star (Baumgarte & Shapiro 1999; Saijo et al. 2002; Carr et al. 1984; Nakamura et 
al. 2001). The seed BH will grow as long as there is a gas supply for the inner disk. This gas 
reservoir is replenished by viscous inflow from the outer disk into the inner region of the disk 
and augmented by the increase in the inner radius due to the growth in the BH mass. A 
sufficient cold gas with transferring rate Mext could be supplied from the quasar host galaxy, 
for instance, mass loss from stars, stellar disruption and collisions etc., if the seed BH is 
embedded in densely populated stellar systems (Tsuruta 1977); and from the progenitor 
galaxies if the seed BH formed in major mergers (Kauffmann et al. 2000). According to the 
model (i) of Cattaneo et al. (1999) that M^xt is proportional to the mass that has formed 
stars in the host galaxy. We take M^xt as a model parameter in this paper. If there is lack of 
any limiting factor (Burkert & Silk 2001), all the gas supply will eventually be accreted and 
fall into the BH. Within the framework of our model, Eddington limit (Small & Blandford 
1992) is used for the growth of BH, and the ionization instability of the quasar disk provides 
the self-regulation that limits the growth of the central BH in three different phases. 


2.1. The S-shaped ionization instability of the quasar disk 

The most realistic scenario for an accretion disk around a BH is time dependent. The 
stationary local disk solutions form a characteristic S-shape on the M — S (accretion rate vs. 
surface density) plane. The accretion disk is thermally and viscously unstable when ^ < 0 
(Meyer et al. 1984). This is the case for accretion rates in the range (Mi, M 2 ), where the 
S-shaped has a negative slope ( Fig.l). A characteristic S-shaped instability consists of three 
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branches (lower, middle, and upper branches) creating three possible accretion phases of the 
quasar disk. In the case of Keplerian motion the surface density evolution may be written 
(Pringle 1981; Smak 1984) as 


(9S 

'm 


3^ 

r dr 

Tir dr 


. 1/2 


dr 

.1/2 


(z/Sr^^) 


+ 


l/2^ 
r ' , 
out} 


dM 

dr 


1 dM 
27rr dr 


( 1 ) 


where rout is the outer radius of the disk and v is the kinematic viscosity coefficient. The 
local accretion rate is given by 

M = 67rr^/^.;^(r^/^z/S) . (2) 

dr 

Equation (1) describes a disk which is continuously supplied not only with the matter but 
also with angular momentum. If the mass inflow is variable with time, the thermally stable 
upper branch of the S-shaped instability (’’active” phase) is characterized by an accretion 
rate higher than the stationary value, while the stable lower branch (’’quiescence” phase) has 
a lower accretion rate than the stationary accretion rate, and so gas is accumulated in the 
outer region of the disk (Siemiginowska et al. 1996). In this model, the disk is assumed to 
be in the conhguration of advection-dominated inflow-outflow solutions (ADIOS)(Blandford 
et al. 1999) in the lower branch, in which most of the material transferred from the quasar 
host galaxy to the accretion disk is ejected in the form of outflowing wind. This is based 
on the observations that the quasar is not active on the lower branch, which hints that 
the efficiency in converting the gravitational energy to radiation is lower, thus the material 
transferred from the host galaxy can be accumulated and clumped into the outer region of the 
disk initially. When the clumped mass in the outer region is transferred into the inner region 
by viscosity, the central seed BH gains mass and begins to grow. At the same time, when 
the surface density of the disk increases to a limit critical value (Figl), the ionization 
instability occurs, the disk jumps from the lower branch to the upper branch. In the upper 
branch, the gas in the disk completely ionized and the disk restores to stable. Unlike the 
lower branch, in the upper branch, the disk is highly active, radiation from the disk reaches 
the highest possible luminosity, restricted only by the Eddington limit. At this stage, the 
observed galaxy is naturally in a highly active phase and the quasar is like a bright quasars 
(Ho 1996). In order to obtain a limit cycle behavior of ionization instability of the disk, 
some stabilizing mechanism must be taken into account in the upper branch. Janiuk et al. 
(2002) proposed that a hot outflowing wind as a result of the disk instability should develop 
to provide a stabilizing mechanism in the upper branch. Energetically, the request for such 
an outflow is consistent with observations (Fender et al. 1999; Reeves et al. 2001). Unlike 
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the outflows of lower branch, the hot outflowing wind of the upper branch have different 
properties, which is probably produced by another mechanism, and the mass-loss rate in the 
upper branch is negligible in comparing with the accretion rate into the BH (Janiuk et ah 
2002). After an outburst, the outer disk attempts to refill the central disk via steady inflow 
and begin the next cycle. The whole cycle repeats roughly on the viscous time of the outer 
disk. Because the recurrence time-scale of the disk is longer, the individual quasar changes 
too slowly to be observed in the Galaxy. It is reasonable to assume that every quasar, as 
well as the growth of the BH, has already experienced a sufficiently large number of cycles 
between the quiescent and active phase in the cosmic time. 


2.2. The available gas mass accumulated into the the quasar disk 


In subsection 2.1, we analyzed the evolution of the quasar disk along S-shaped instability 
and how the BH gains its mass. It is interesting to calculate quantitatively the BH growth 
on three phases in one cycle. We first estimate that how much mass from the quasar host 
galaxy can be transferred and accumulated into the outer region of the disk. Considering our 
assumption that the disk in the lower branch is in the ADIOS configuration and the mass-loss 
of the wind from the quasar disk in the upper branch can be negligible (Janiuk et ah 2002), 
then the transferred mass from the host galaxy to the disk thus can be divided into two 
components when quasar disk evolves in one cycle: one component is the infalling material 
with mass {AMbh)acc, which is also the available gas mass accumulated into the outer region 
of the quasar disk around the BH ; another is outflowing material with mass AM^ut which 
is ejected into the ISM in which quasar lies. To calculate AMout and {AMhh)acc, we use the 
pseudo-Newtonian RADIOS model proposed by Becker et al. (2001) to describe the relation 
between the energy emitted into the wind and the amount of mass per unit radius per unit 
time, which is 
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out 
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where e must satisfy max < e < Dq, which ranges 0.025 < e < 0.1 with Dg = 

0.1, Do is positive constant, r is the radius expressed in Schwarschild radii [R^ = ^^^^), 
Ibh = ^ , Lbh is the bolometric luminosity emitted by material accreting onto 

the BH, and Mout is the mass-loss rate. MEdd and LEdd is the Eddington accretion rate and 
luminosity, respectively. Tin and r-o^t are the inner and outer radius of the disk, respectively. 
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For this reference model, we adopt = 6 and rout 10^ (Narayan et al. 1996, 1997a), 
which gives the numerical value of integral in Eq (3) ~ 0.03. 

In the case of ADIOS situation, the accumulated mass in the outer region of the quasar 
disk and the outflowing mass on viscous time-scale {tyis) in one cycle can be estimated as 

(AMbh) acc acctvis (4) 

AMout ~ Moutt vis •) (5) 

where {Mbh)acc is the accumulated mass accretion rate in the quasar disk around the BH. 
Eventually, the accumulated material on the quasar disk should fall into the BH. Taking 
Lbh ~ 0.1{Mbh)accC^, and combining Eq.(3) and Eq.(4), we obtain 

{AMbh) acc ^ 0.03eiifoMtf vis • (6) 

Eq.(6) gives the available maximum mass in the outer region of the disk in one cycle. Let 
the values of e varies as 0.025 < e < 0.1, we get 

7.5 X IQ-^MoutUis < {AMbh)acc < 3 X IQ-^Mouttuis- 

We conclude that {AMbh)acc is about 10“^ of the total outflowing mass. Combining Moxtivis = 
AMout + {AMbh)acc and Eq.(5), we thus have 

AMout ~ Moxtt vis 1 (7) 

and then 

{Mbh)acc ~ O.OSeMext ■ (8) 

We shall treat M^xt as a model parameter, which is determined by the total stars formed in 
the host galaxy (Cattaneo et al. 1999). 


2.3. The mass growth of the black hole along the S-shaped instability in one 

cycle 

In this section, we study the growth of the central BH along the S-shaped instability 
in one cycle. We argue that in order to realize the next cycle of S-shaped instability, the 
accumulated mass {AMbh)acc of the quasar disk must be either swallowed by the BH or 
ejected away in the upper branch when the quasar is in highly active state. As the upper 
branch of the S-shaped instability involves physics that is currently opened (Janiuk et al. 



2002; Blandford 2002; Mirabel et al. 1999). We estimate the properties of the outflowing and 
infalling in the upper branch according to the observations. It is known that jets are among 
the most ubiquitous of astrophysical phenomena, and are associated with objects ranging 
from proto-stars and binary systems in our Galaxy to SMBHs in the centers of AGNs. Large 
asymptotic Lorentz factors are implied by observations of superluminal motion in many 
blazars (Dermer et al. 1992; Vermeulen et al. 1994), and the EGRET instrument on board 
the Gompton Gamma Ray Observatory (GGRO) has also detected intense y-ray flares from 
many of these sources (Wehrle et al. 1998; Kataoka et al. 1999). The mechanisms responsible 
for producing the observed jets are still poorly understood, bnt collimated ontflows in general 
appear to be associated with objects that derive their Inminosity from the accretion of matter 
onto a gravitating object (Dermer et al. 1992). The observations of relativistic jets shows 
that asymptotic Lorentz factors is about 'jjet > 10 (Meszaros, etal. 1992; Quinn et al. 1996). 
For an accretion rate Mfe/j, the emitted luminosity can be estimated as 

Lbh = O.lMfohC^ • (9) 


If the ontflow gas expands adiabatically, the total energy per unit mass in the jet at a rate 
Mjet follows 

Ljet = IjetMjetf? , . ( 10 ) 


It therefore follows immediately that 


Mjet 

Mbh 


1 


< 1 . 


( 11 ) 


Eq.(ll) thus shows that the mass flux in the jet on the upper branch is only less than 3% 
of the mass flux into the BH. This simple energetic argument predicts that ontflowing mass 
can be ignored in comparing with inflowing mass in the npper branch. We therefore neglect 
the mass loss in the npper branch. This means that all accumulated mass onto the quasar 
disk will be swallowed by the BH before a new cycle starts. Gonseqnently, we have 


= {AMbh)acc = O.OSeMexttvis , 


( 12 ) 


where AMbh is the total mass growth of the BH gained in each cycle of the S-shaped insta¬ 
bility. 

We are also interested in the detail, with observation constraints, how the BH grows up 
along the S-shaped instability of the quasar disk. It is known that the surface density and 
the local accretion rate of the quasar disk evolves on the viscous timescale t^s (Lightman 
1974) tvis ~ R/vr (where vr is the radial velocity, and R is the accretion radius of the disk). 
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thus the life time of a quasar accretion disk in each cycle is tms- The recurrence timescale 
^current of S-shaped instability of the quasar disk is addressed as tcurrent ~ t^s, which operates 
over 10*^ - 10® yr (Siemiginowska et ah 1996). Time-scales spent by the quasar disk in lower, 
transition and upper branches in one cycle are tiow, Uran and tup, respectively, which are 
given by Siemiginowska et ah (1996), 


^low 


( 13 ) 

tup 

~ 0.\tyis , 

( 14 ) 

'tran 


( 15 ) 


We assume that the total growth mass of the BH in one cycle is built up when 

the quasar disk passes three different phases. Then the total growth mass of the BH in one 
cycle is 

AMi)fi M^hitiou) T Mhfiutup T , (16) 

where Mbhi, M^ht and Mbhu are the accretion rate of the BH in lower, transition and upper 
branches, respectively. Each accretion rate is limited by Eddington accretion rate. 

The accretion rate in the transition phase can be estimated from the feature of S- 
shaped instability. It is noted that the ionization instability of the quasar disk occurs when 
the accretion rates of the disk are in the range (M 2 , Mi) (see, Fig.l). We thus approximate 

Mbht = - (^Mbhi + Mbhu^ , ( 17 ) 

Inserting Eq.(13) - Eq.(15) and Eq.(17) into Eq.(16), we get 

AMbh ~ 0.75Mbhltvis + OAMbhutvis 

-1-0.075 (^Mbhi + Mbhv^ tuis ■ (18) 

Comparing Eq.(12) with Eq.(18), we get 

0.03eMe,i = t).S2bMbhi + 0.175Mfe;,„ . (19) 

If 0.01 < <0.1 (Mirabel et ah 1999), the accretion rates in three phase can be estimated 

roughly 

1.64 X 10“^eMe„i < Mbhi < 1-16 x 10“^eMea,i 

1.16 X 10-^eMe,i < Mbhu < 1.64 x 10-^eMe,t 

2.28 X lO-^eMexi < Mbht < 6.38 x lO'^eMe^t 


( 20 ) 
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If Mbhu = MecIcI, the limits of Mext required by our model can be inferred from Eq.(20), 
which gives 

<Mext< 86 MqJt~^ (21) 

for e = 0.1 , 

24.4Moyr-i < Mext < 345 MqJt-^ (22) 

for e = 0.025 . 

The mass growth of the BH in three phases of one cycle can be obtained 

1.23 X 10-\Mexttvis < ^Mbhl < 8.7 X IQ-hMexttvis 

1.16 X IQ-hMexttvis < ^Mbhu < 1.64 X IQ-hMexttvis 

3.42 X 10~^eMexttvis < ^Mhht < 9.57 x eMextUis 

(23) 

where AM^hi = MbhiUow, AM^ku = Mbhutup, and AMbht = Mbhtttran are the mass growth of 
the BH in the lower, upper and transition branches, respectively. 


3. The present mass of BH and the BH - bulge mass relation 

The total mass of the resulting BH at present is the sum of the growth mass through 
the seed BH accreting in each cycle during the Hubble time. Taking into account model 
assumption that the accretion rate onto the quasar disk is limited by Eddington-limits, we 
can rewrite Eq.(12) 

dMbh _ f 0.03eMexi ; if O.OSeMext < MEdd 
dt \ MEdd ; if O.OSeMext > MEdd 

The solutions of Eq.(24) are: 

(i) If O.OSeMext A MEdd, 

Mhhif) = O.OSeMextt + Mbhiff), (25) 

where t ranges 0 < t < tn, tu is the Hubble timescale (tn = lOGyr). Mbhif) and Mbh{0) 
are the present BH mass and the seed BH mass, respectively. 

For 0.025 < e < 0.1, the seed BH mass Mbh{0) are: 


2 X 10®Mo < Mbh{0) < 3 X IO^Mq 


(26) 
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Then the present BH mass becomes 2 x 10 ®Mq and 3 x lO^M© for the seed BH mass 2 x IO^Mq 
and 3 x IO^Mq respectively. 

(ii) If 0.03eMea;t > MEdd , 

M,n{t) = , (27) 

where tEdd is Eddington accretion time. This solution is not important because the BH 
will Erst grow exponentially, quickly O.OSeMe^t < MEdd{Mhh(t)) it switches back to the first 
solution ( Eq.(25)). Since the time spends in the exponential growth is less than the Hubble 
time so this stage can be ignored unless the seed BH mass is extremely small. Consequently, 
the solutions of Eq.(24) suggest that a typical quasar BH with a final mass 2 x 10 ®Mq can 
be grown from a seed BH with mass 2 x IO^Mq. 

As assumed that the material transferred from the host galaxy is proportional to the 
mass that has formed stars, this situation is strictly true in model (i) of Cattaneo et al. 
(1999), then AMout ~ Mexttvis oc M^uige, where Mbuige is the mass of host galactic bulge. 
Taking into account of Eq.(7) and Eq.(12), and let 0.025 < e < 0.1, we thus obtain the mass 
ratio of the BH to the galactic bulge 

7.5 X 10"^ < < 3 X 10“^ . (28) 

^bulge 


4. Discussion and conclusion 

A quasar disk model associated with the ionization instability is proposed to account 
for the build up of SMBHs in quasar host galaxies. The growth of the BH is studied when 
the quasar disk evolves along S-shaped instability. We shows that the main growth of 
the BH occurs in the upper branch when quasar disk is in a highly active state. Two 
model parameters e and Mext are involved. Though in quasar accretion disks the gas supply 
mechanism is unknown (Siemiginowska et al. 1996), the limits for an available gas supply 
with 6 MQyr“^ < Mext A 86 M 0 yr“^ from the quasar host galaxy is deduced based on 
observations and the limits of Eddington-accretion rate. Our conclusion can be summarized 
as: (i) We first demonstrate that the BHs of the quasar host galaxies can grow up through 
accretion with two model parameters. With 0.025 < e < 0.1, a BH with mass 2 x lO^M© can 
grow up from a seed BH with mass 2 x IO^Mq. (ii) Secondly, the observed BH - bulge mass 
relation (Magorrian et al. 1998) can be inferred, which is 7.5 x 10“^ < < 3 x 10“^ in 

^'^■^bulge 

our model parameters, (hi) Furthermore, in present scheme, the massive seed black hole is 
required in order to build up super-massive BH through accretion if the gas is provided from 
the quasar host galaxy. This hints that a low-mass stellar BH may not grow to super-massive 
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black hole unless without Eddington-limit. The details for the formation of massive seed BH 
(~ IO^Mq) and if there is a sufficient cold gas supplied to the BH are beyond the scope of 
this paper. 
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